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Resonant interaction between two (or more) progressive interfacial waves produces ex- 
ponentially growing instabilities in idealized, homogeneous or density stratified, inviscid 



shear layers (Holmboe 1962 1 . Resonance occurs when the two waves attain a "phase- 



locked" configuration. In this paper we have generalized the mechanistic picture of shear 



instabilities described in Holmboe (1962). Unlike Holmboe, we do not initially assume 
the wave type (e.g. vorticity wave or gravity wave), nor do we impose the normal- mode 
waveform (which only accounts for exponential growth) . Starting from first principles, we 
demonstrate that two oppositely propagating interfacial waves, having arbitrary initial 
amplitudes (which are small enough to satisfy linearity) and phases, eventually phase- 
lock, provided they satisfy a certain condition. We show it to be the necessary and 
sufficient condition for exponentially growing instabilities in idealized shear layers. We 
investigate three types of shear instabilities - Kelvin Hclmholtz, Holmboe and Taylor. It 
is shown that the above-mentioned condition provides a range of unstable wavenumbers 
for each instability type, and this range matches the predictions of the canonical normal- 
mode based linear stability theory. The wave interaction problem is interpreted using 
dynamical systems approach. It is shown that the dynamical system is of "source-sink" 
type, and its equilibrium points represent the growing and decaying normal-modes of the 
discrete spectrum. Once arbitrarily initialized, the two waves adjust each other's ampli- 
tudes and phases until they reach the resonant configuration, which corresponds to the 
growing normal-mode of the spectrum. The pre-resonance stage characterizes the non- 
modal instability of the shear layer, favoring rapid transient growth processes. Optimal 
growth is found to occur when the the two waves are in quadrature. 



1. Introduction 

Statically stable density stratified shear layers are ubiquitous in the atmosphere and 
oceans. Such shear layers can become hydrodynamically unstable, resulting in turbulence 
and mixing in geophysical flows. Turbulence and mixing strongly influence the atmo- 
spheric and oceanic circulation - processes known to play key roles in shaping the weather 
and the climate. Hydrodynamic instability, the precursor of turbulence, is characterized 
by the growth of wavelike perturbations in a laminar base flow. Such perturbations can 
grow at an exponential rate, transforming the base flow from a laminar to turbulent state. 
In the present study, we will investigate theoretically the underlying mechanism(s) lead- 
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ing to the exponential growth of small wavelike perturbations in idealized homogeneous 
and stratified shear flows. 

The classical method used to determine flow stability is the normal-mode approach 
of linear stability analysis ( Drazin fc Reid|2004 ). Under the normal-mode assumption, a 
waveform can grow or decay but cannot deform. Normal-mode perturbations are added 
to the base flow, followed by linearizing the governing Navier Stokes equations about the 
laminar background state. For inviscid, density stratified shear flows, the normal-mode 
formalism leads to the Taylor- Goldstein equation, derived independently by Taylor ( 1931 1 
and Goldstein (1931). The Taylor- Goldstein equation is an eigenvalue problem which 
calculates the wave properties like growth rate, phase-speed, and cigenfunction associated 
with each normal-mode. For stability analysis, the range of unstable wavenumbers and 
the wavenumber corresponding to the fastest growing normal-mode are of prime interest. 
Generally, the fastest growing mode dominates over all other modes, therefore it governs 
the behavior of the transitional flow. In practice, the Taylor- Goldstein equation is found 
to accurately capture the onset of instability ( Thorpe||1973 ), and it also provides a first 



order description of the developing flow structures (Thorpe 1973 Tedford et «/.| |2009). 



The normal-mode approach to linear stability analysis has a few principal shortcom- 
ings. Firstly, it is an involved mathematical exercise providing little insight into the 
physical mechanism(s) responsible for hydrodynamic instability. The answer to why an 
infinitesimal perturbation vigorously grows in an otherwise stable background flow is pro- 
vided in the form of non-intuitive mathematical theorems - Rayleigh's inflection point 
theorem and Fj0rtoft's extension for the case of homogeneous flows, and Miles-Howard 
criterion for stratified flows ( Drazin fc Reid|2004[ ). Since linear instability is the first step 
towards understanding the more complicated and highly elusive non-linear processes like 
chaos and turbulence, it is desirable to formulate alternative theories which are able to 
provide intuitive explanations. A second drawback of the normal-mode approach is the 



normal-mode assumption itself. The extensive work by Farrell (1984), Trefethen et al. 
( 1993 ), Schmid & Henningson ( 2001[ ) and others have shown that shear allows rapid non- 
modal transient growth due to non-orthogonal interaction between the modes. |Farrell fe"| 



Ioannou (1996) developed the Generalized stability theory for linear dynamical systems, 



and showed the process of obtaining the optimal non-modal growth from a singular value 
decomposition of the propagator matrix of the linear dynamical system. 

Theories have been proposed to provide an intuitive understanding of the hydrody- 
namic instability process. Probably the first mechanistic picture of stratified shear insta- 
bilities was provided by Holmboe (1962). Using idealized velocity and density profiles, 
Holmboe postulated that the resonant interaction between stable propagating waves, 
each existing at a discontinuity in the background flow profile (density profile discontinu- 
ity produces gravity waves and vorticity profile discontinuity produces vorticity waves), 
yields exponentially growing instabilities. He was able to show that Kelvin Hclmholtz 
instability (Rayleigh 1880) is the result of the interaction between two vorticity waves 
(also known as Rayleigh waves). Moreover, Holmboe also found a new type of instability, 
now known as the "Holmboe instability" , produced by the interaction between vorticity 
and gravity waves. Bretherton, a contemporary of Holmboe, proposed a similar theory to 
explain mid-latitude cyclogenesis ( Bretherton|1966 ) . He hypothesized that cyclones form 
due to a baroclinic instability caused by the interaction between two Rossby edge waves 
(vorticity waves in a rotating frame of reference), one existing at the earth's surface and 
the other located at the atmospheric tropopause. The theory proposed by Holmboe and 



Bretherton has been refined and re-interpreted over the years, see Cairns ( 1979 1; Hoskins 



et al. 


(1985 


); 


Caulfield 


(1994 


); 


ter et al. 


( 


2012 


). As reviewed 
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two edge waves in an idealized homogeneous or stratified shear layer occurs when these 
waves attain a phase-locked state, i.e. they are at rest relative to each other. Maintaining 
this phase-locked configuration, the waves grow equally at an exponential rate. 

A very different mechanistic picture for describing shear instabilities was proposed 
by Lindzen and co-authors (summarized in the review article by |Lindzen| ([1988)). This 
theory, known as the "Over-reflection theory" , proposes that under the right flow con- 
figuration, over-reflection of waves can continuously energize an advective "Orr process" 
( |Orr|190"7| ) which is finally responsible for the perturbation growth. In this paper we will 
only consider the wave interaction based interpretation of shear instabilities. In our opin- 
ion, the wave interaction mechanism provides a more intuitive understanding because of 
its kinematic approach. 



In recent years, Heifetz and co-authors (Heifetz et al. 1999 2004 Heifetz & Methven 



2005 ) have extensively studied the interaction between Rossby edge waves. By not lim- 



iting the Rossby edge waves to be of the normal-mode type, they were able to obtain 
non-modal instability and transient growth mechanisms. While Heifetz et al. (2004) de- 
rived the governing equations using the Hamiltonian approach, Heifetz & Methven (2005 1 
used the stream function- vorticity approach to do the same. Their successful attempt has 
motivated us to formulate a generalized interfacial wave interaction technique for study- 
ing homogeneous and stratified shear instabilities. Without forcing the wave to be of the 
normal-mode type, and furthermore, without assuming any particular type of waveform 
(e.g. gravity wave or vorticity wave), we have formulated the governing equations for wave 
interaction. Our equations are derived from the linearized kinematic and dynamic (for 
stratified flows) conditions. Unlike [Heifetz et al. (2004) and |Heifetz fc Methven (2005), 
our key variables are the vertical displacement and the vertical velocity of the wave. 
The choice of variables, along with the derivation methodology, provide a deeper and 
probably more intuitive understanding of the wave interaction process. 

The outline of the paper is as follows. In f[2]we provide theoretical background of simple 
progressive waves, and focus on two types of waves - vorticity waves and internal gravity 
waves. The wave theory in this Section is more generalized than usually reported in the 
literature. In £j3] we investigate the mechanism of interaction between two progressive 
waves. We undertake a dynamical systems approach to better understand the wave in- 
teraction problem, especially the resonant condition. In this wave interaction perspective, 
we describe three types of shear layer instabilities in Q - Kelvin Helmholtz instability 
(resulting from the interaction between two vorticity waves), Holmboe instability (re- 
sulting from the interaction between a vorticity wave and an internal gravity wave) and 
Taylor instability (resulting from the interaction between two internal gravity waves). 



2. Linear wave(s) at an interface 

Let us consider an inviscid, incompressible, two dimensional flow with hyd rostatic 
variation of pressure. In such flows, the evolution of vorticity $7 is given by (Turner 



1973): 



(2.1) 



Dn_dQ _dCl „dn_gdp 

Dt dt dx dz p dx 
Here u, w and p respectively imply horizontal (streamwise) velocity, vertical velocity 
and density. The last term of Eq. (2.1 1 signifies the baroclinic generation of vorticity. 
We assume the mean flow to depend only on the vertical coordinate z. Furthermore, we 
decompose the flow field into mean and fluctuating components: VL = D,(z) + Cl(x, z,t), 
u = U(z) + u(x, z, t), w — w(x, z,t), p — p(z) + p(x, z,t). Substituting them in Eq. (2.1 1 
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and ignoring the non-linear terms, we get 



DQ 

T5t 



~dt 



U 



dx 



-w 



dz 



g_dp_ 
p dx 



(2.2) 



Here D/Dt implies the linearized material derivative. Eq. (2.2) reveals that perturba- 
tion vorticity Cl is generated from two distinct sources - vertical gradient of background 
vorticity (dQ/dz) and horizontal gradient of density perturbation (dp/dx). 

In the present study we consider multi-layered flows with constant density and vorticity 
in each layer. This configuration makes the equations of motion for perturbations within a 
layer to be the same as that in an irrotational background flow. The interface between two 
adjacent layers signifies a discontinuity in vorticity or density. The former is a vorticity 
interface, while the latter is a density interface. Let such an interface existing at a location 
z — Zi be perturbed by an infinitesimal vertical displacement rji 1 given as follows: 



Vi = ^(t) cos[ 



<M*)] 



(2.3) 



This displacement manifests itself in the form of stable, progressive wave(s), the ampli- 
tude and phase of which are A Vi and (\> m respectively. While a vorticity interface produces 
a vorticity wave, two oppositely traveling gravity waves are produced in the case of a 
density interface. We have assumed the interfacial displacement (or the wave) to be 
monochromatic, having a wavenumber a. Moreover, the interface satisfies the kinematic 
condition - a particle initially on the interface will remain there forever. The linearized 
kinematic condition is given by 



dt 



dm 

dx 



(2.4) 



where Ui = U(zi) is the background velocity and Wi is the vertical velocity at the 
interface. We prescribe the latter to be as follows: 



Wi = A w . {t) cos [ax + 4> w 

Here A Wi is the amplitude and <f> Wi is the phase of w, 
up a Laplacian velocity field (u,w) away from the interface (Holmboe 1962 Caulfield 
1994), the vertical component of which is given by: 



>(*)] (2-5) 
. The interfacial deformation sets 



-a\z—Zi | 



W; 



(2.6) 



Eq. (2.6) indicates that the vertical velocity disappears at infinity, a condition known as 



wave evanescence 



Substituting Eqs. (2.3) and (2.5) in Eq. (2.4), we obtain 

i {ax + (f> Vi ) — A m (^aUi + 4> rji ^ sin (ax + <fi rii ) = A Wi cos (ax + 4> Wi ) 



A Vi cos ( 



(2.7) 



Here 



G [— 7r, 7r]. The frequency and the growth rate of a wave are respectively 



defined as Oi = —<j> Vi and ji = A Vi /A Vi (overdot denotes d/dt). Using these definitions 
in Eq. (2.7), we get 



<7; = alii - Wi sin (A<f>a) 
7i = Wi cos (A&j) 



(2.8) 
(2.9) 



where A<fiu = <p Wi — (p m , Eq. (2.8) shows that the frequency of a wave consists of two 
components - (i) the Doppler shift alJi due to the background velocity, and (ii) the 
intrinsic frequency — C0j sin (A^), where Wi = A w ./A Vi . The phase-speed Cj = <7j/a of 
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the wave is found to be 



= Ui- 



sin (A<j>u) 



(2.10) 



The last term (including the negative sign) denotes the intrinsic phase-speed.Notmg that 
a wave in isolation cannot grow or decay on its own, Eq. (2.9) demands that |A^jj| = 
7r/2. Therefore for a stable wave, the vertical velocity field at the interface has to be 



in quadrature with the interfacial deformation. According to Eqs. (2.8) and (2.10), the 



quadrature condition makes the magnitude of the intrinsic frequency and the intrinsic 
phase-speed to be Wi and Wi/a respectively. The intrinsic direction of motion of the 
wave, however, is determined by A<^. For waves moving to the left relative to the 
interfacial velocity Ui, Ac6^ = tt/2. Similarly for right moving waves, A(j>u — — tt/2. 
When such a stable, progressive wave is acted upon by external influence(s) (e.g. when 
another wave interacts with the given wave, as detailed in |3j, the quadrature condition 
is no longer satisfied, i.e. |A0ii| ^ tt/2. Therefore, the wave may grow (7, > 0) or decay 
(~fi < 0), and its intrinsic frequency and phase-speed may change. 

In our analyses we will consider two types of progressive interfacial waves - vorticity 
waves and internal gravity waves. 

2.1. Vorticity Waves 

Vorticity waves, also known as Rayleigh waves, exist at a vorticity interface (i.e. regions 
involving a sharp change in vorticity). Such interfaces are a common feature in the 
atmosphere and oceans. In a rotating frame, the analogue of the vorticity wave is the 
Rossby edge wave which exists at a sharp transition in the potential vorticity. When 
Rossby edge waves propagate in a direction opposite to the background flow, they are 
called "counter-propagating Rossby waves" or CRWs ( |Heifetz et a/TpM)] ). 

In order to evaluate the frequency ix>i of vorticity waves, let us consider a velocity 
profile having the form 



U(z) 




(2.11) 



Here the constant S = Ui/zi is the vorticity, or the shear in the region z ^ z^ (Carpenter 
et al. 2012|. Eq. (2.11) shows that the vorticity dU/dz is discontinuous at z = Zi. This 
condition supports a vorticity wave. A deformation rji of the interface adds vorticity S 
to the upper layer and removes it from the lower layer, creating a vorticity imbalance 
and thereby a mechanism for wave propagation. The horizontal component Ui of the 
perturbation velocity field set up by the interfacial deformation undergoes a jump at 
the interface, the value of which can be determined from Stokes' Theorem (Eq. (3.2) of 



Holmboe (1962)): 



Sf]i 



(2.12) 



By taking an x derivative of Eq. (2.12) and invoking the continuity relation, we get 

dw~ 



dz 



dz 



= S 



drjj 

dx 



By substituting Eqs. (2.3) and (2.5) in Eq. (2.13), we obtain 

S 



w, 



S sin (ax + ^ ) 



2 cos (ax + <j) 
The fact that a>i is always positive demands 



2sin(A<fe) 



(2.13) 



(2.14) 



Ac 



-sgn(S) 



(2.15) 
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where sgn( ) is the sign function. From Eq. (2.14|, the intrinsic frequency of a vorticity 
wave is found to be —5/2. The phase-speed Cj can be evaluated by substituting Eq. ( 2.14 1 
into Eq. \2.1Q\ : 

c t = Ut - ~ (2.16) 
2a 

If 5 > 0, the vorticity wave moves to the left relative to the background flow. Alternative 
derivations of the frequency and phase-speed of a vorticity wave can be found in several 



Caulficld ( 


1994), 


Sutherland ( 


20101 or 



2.2. Interfacial Internal Gravity Waves 

Intcrfacial gravity waves exist at a density interface, i.e. regions involving sharp change 
in density. The most common example of interfacial gravity wave is the surface wave 
existing at the interface of air and water. Here we will be considering interfacial internal 
gravity waves (hereafter, gravity waves) only. Such waves exist in density stratified flows 
having a thin density interface (pycnocline) . Since most natural water bodies like lakes, 
estuaries and oceans are density stratified, gravity waves are virtually ubiquitous. 

In the case of gravity waves, the intrinsic frequency — tUjSin (A(fra) can be evaluated 
by considering the dynamic condition. The latter implies that the pressure at the density 
interface must be continuous. Let the density of upper and lower fluids be p± and pi 
respectively. Then the linearized dynamic condition at the interface after some simplifi- 



cation becomes (Eq. (3.13) of Cauffield (1994)) : 



dt 



U, 



dx 



9 drjj 
2a dx 



(2.17) 



Here g = g(p2 — Pi)/ Pa is the reduced gravity and po is the reference density. Under 
Boussinesq approximation po ~ pi ~ pi- By taking an x derivative of Eq. (2.17) and 
using the stream function relation {ui,Wi} 



dwj 
~dt 



+ U, 



{—dipi/dz, dipi/dx}, we get 

dwj _ g d 2 ^ 
dx 2a dx 2 



(2.18) 



Substitution of Eqs. and ([2^ in Eq. \2.18\ yields 

A Wi cos (ax + (j) Wi ) — A Wi <p Wi sin (ax + <f> Wi ) 

g a 



-aUiA w . sin (ax + 4> Wi ) 



-A Vi cos (ax + 0,,J 



(2.19) 



The quantity 
obtain 



w 



aC{. On substituting this relation in Eq. (2.19), we 

g sin (A(j)ii) 



2 (Ui-ct) 

Since o>j is a positive quantity, Eq. (2.20) demands 



A<fe = -sgn(Z7i - a) 



(2.20) 



(2.21) 



An important aspect of Eq. (2.20) is that it has been derived independent of the 



kinematic condition. The presence of single or multiple interfaces does not alter the 



expression in Eq. (2.20), implying that this equation provides a generalized description 
of cu;. Inclusion of kinematic condition yields an expression for a>i which is simpler but 
problem specific. For example, when a single interface is present, inclusion of kinematic 



condition in Eq. (2.201 produces the well known expression for gravity wave frequency. 
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Substituting Eq. (2.8) (this equation has been derived from the kinematic condition for 



a single interface, i.e. Eq. ( 2.4 1 ) in Eq. (2.20) and considering only the positive value, we 



obtain the dispersion relation for gravity waves: 




(2.22) 



Moreover Eq. (2.9) requires |A<^j| = w/2. Substituting it along with Eq. (2.20) in Eq. 
(2.10)) produces the well known expression for the phase-speed of a gravity wave: 



- Ui± 



(2.23) 



The above equation shows that each density interface supports two gravity waves, one 
moving to the left and the other to the right relative to the background velocity Z7j. 
Alternative approaches to deriving the frequency and phase-speed of a gravity wave can 



be found in Caulfield (19941; Sutherland (20101; Carpenter et al. (2012 1). 



3. Interaction between two linear interfacial waves 

Let us now consider a system with two interfaces, one at z = Z\ and the other one at 
z — Z2- The linearized kinematic condition at each of these interfaces is then given by: 



dm v drg 

dt dx 

dm LT dm 

dt dx 



e -«\*l~'2\ V]l+W2 



(3.1) 
(3.2) 



It has been implicitly assumed that both waves have the same wavenumber. The R.H.S. 



of Eqs. (3.1 )-(3.2 1 reveal the subtle effect of wave interaction, and can be understood as 



follows. The effect of u>i extends away from the interface z\, hence it can be felt by a 
wave existing at another location, say z 2 . Therefore the vertical velocity of the wave at z^ 
gets modified - it becomes the linear superposition of its own vertical velocity W2 and the 
component of w\ existing at Z2 ■ This phenomenon is also known as "action-at-a-distance" , 
see 



Heifetz & Methven (2005) 



On substituting Eqs. (2.3 1 and (2.5 1 in Eqs. (3.1)-(3.2), we get 



A Vl cos (ax + 4> m ) — A, n [aUi + <j> ni \ sin (ax + (f> m ) — 
A Wl cos (ax + (j> Wl ) + e~ a \ Zl ~ Z2 \A W2 cos (ax + (f> W2 ) 

A V2 cos (ax + (f> V2 ) — A V2 (aU2 + sin (ax + 4>., l2 ) = 
e~ a ' zl ~ Z2 'A wl cos (ax + <j> Wl ) + A W2 cos (ax + <f> W2 ) 



(3.3) 



(3.4) 



Proceeding in a manner similar to |2l the growth rate 7* and phase-speed c, of each wave 
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is found to be 



7i 



A, : 



■ cos (A0 n ) + ^l e - Q l Zl - Z2 l cos (A&a) 



ci = C/i 

■An 



4^ sin (A0 n ) + ^ e -"^-^ sin (A0 12 ) 



72 



COS (A0 22 ) + ^l e - Q l Zl - Z2 l COS (A02!) 



c 2 = U 2 

a 



^V2 

^ sin (Afe) + ^?L e -«l»i-».l sin (A0 21 ) 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



Here A0jj — When — 22] — > oo, the two waves get uncoupled, and we 

recover Eqs. (2.9)- (2.10) for each wave. As argued in Sj2j a wave in isolation cannot grow 
or decay on its own. Therefore, the first term in each of Eq. (3.5) and Eq. (3.7) should 
be equal to zero, implying |A<^j| — tt/2. 

In all our analyses, we will be considering a system with a left moving top wave ( A</>n = 
7r/2) and a right moving bottom wave (A<^ 22 = — 7r/2), the wave motion being relative to 
the background velocity at the corresponding interface. Let the phase shift between the 
bottom and top waves be $ = (f> V2 — (f> Vl . Therefore A</>i 2 = $ — tt/2 and A</> 2 i = tt/2 — <f>. 
Defining amplitude ratio R = A m /A ri2 , we re-write Eqs. (3.5)-(3. 



to obtain 



71 = ^e- a l Zl - Z2 lsin$ 
R 



cx = Ux- 



1 



a 



w 1 -^-e-^ Zl - z ^cos^> 
R 



7 2 = i?cuie" Q|2l ~ 221 sin$ 

c 2 = U 2 + - 
a 



a> 2 - i?a)ie" a|21 " 221 cos$ 



(3.9) 

(3.10) 
(3.11) 
(3.12) 



Eqs. (3.9)-(3.12) describe the linear hydrodynamic stability of the system. Unlike the 



conventional linear stability analysis, we did not impose normal-mode type perturbations 
(they only account for exponentially growing instabilities) in our derivation. Therefore 
the equation set provides a non-modal description of hydrodynamic stability in idealized 
(multi-layered) shear flows. We refer to this theory as the "Wave Interaction Theory 
(WIT)". WIT is only applicable to those hydrodynamic stability problems where the 
discrete spectrum dynamics is of interest and the continuous spectrum can be neglected. 
A schematic description of the process of wave interaction is illustrated in Fig. [I] 



Subtracting Eq. (3.11) from Eq. (3.9) and Eq. (3.12) from Eq. (3.101, we find 



— = R ( 7l - 72 ) = (cu 2 - R 2 Wl ) e - a ^- z ^ sin$ 
at 



dt 



= a (ci - c 2 ) = a (Ui — U 2 ) - 



Rw-\ 



cu 2 



S -"l*l-*2| cos$ 



(3.13) 
(3.14) 



The two parameters have the following range of values: R £ (0, oo) and $ 6 [— n, tt]. Eqs. 
(3.13 1- (3. 14 ) represent a two dimensional, autonomous, non-linear dynamical system. The 
two equilibrium points of the system, found by imposing the conditions dR/dt — in 
Eq. ( |3T3[ ) and d$/dt = in Eq. ( |3.14[ ), are given by 



(3.15) 



Wave Interaction Theory 



9 




Z2- 









1 

















Figure 1: Schematic of the interfacial wave interaction mechanism. The deformation 
and associated vertical velocity of each wave is shown by the same color. Interaction 
imposes an additional vertical velocity (shown by different color). The horizontal arrow 
associated with a wave indicates the intrinsic wave propagation direction. Both the waves 
are counter-propagating (move against the background velocity at that location). 



where 



Rum. 




± cos 1 



CUi + CU 2 -a(^l-^2) \ e a\zi-z 2 \ 



(3.16) 
(3.17) 

, (3.18) 

The linear behavior of the dynamical system around the equilibrium points is of inter- 
est. To understand this behavior, we evaluate the Jacobian matrix, J at the equilibrium 
points: 

sin(±0 nm ) 



Eq. (3.17) reveals that the equilibrium points exist only if 

>i + u) 2 -a(Ui- t/gH a \ Zl -z 2 \ 



< 1 



J (Rnm, ±8nm) = -2y/WiW 2 e 



-a\z\— z 2 \ 







sin (±9 nm ) 



(3.19) 



Eq. (3.191 shows that the two eigenvalues corresponding to each equilibrium point are 
equal. Further analysis reveals that every vector at the equilibrium point is an eigenvector. 
The equilibrium point {R n rmS nm ) produces negative eigenvalues, while the eigenvalues 
corresponding to (i? nm , —0 nm ) are positive. Hence the dynamical system represented 
by Eqs. ( 3.13 )-( 3.14) is of "source-sink" type. In terms of the classical normal-mode 



analysis, each equilibrium point corresponds to a normal-mode of the discrete spectrum - 
{Rnm,Qnm) corresponds to the growing normal-mode (signifying exponential growth) and 
{Rnrm—Snm) corresponds to the decaying normal- mode (signifying exponential decay). 



Normal-mode type instabilities can exist only if the condition in Eq. (3.18) is satisfied 



Therefore Eq. (3.18) denotes the necessary and sufficient condition for exponentially 
growing instabilities in idealized, homogeneous and stratified shear layers. 

WIT allows understanding hydrodynamic instability from two different perspectives - 
wave interaction and dynamical systems. According to the former, exponentially growing 
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(a) 




Figure 2: Any initial condition (i?o, ( f ) o) finally yields the resonant configuration 
(Rnm,&nm), provided Eq. (3.18 1 is satisfied. The case depicted here is Kelvin-Helmholtz 
instability (interaction between two vorticity waves) corresponding to a = 0.4. Any other 
shear instability will show qualitatively similar characteristics, (a) $ versus t correspond- 
ing to <!>o = — TT, — 7r/2, 0, 9 nm , tt/2 and tt. The value of Rq is held constant, and is equal 
to 2. (b) R versus t corresponding to Ro — 1/2, l(R nm ), 2 and 5. The value of $o is held 
constant, and is equal to —tt/2. 



instabilities signify resonant interaction between the two waves. From dynamical systems 
point of view, resonance implies "equilibrium condition" (d^/dt = and dR/dt = 0). 
Wave interaction interpretation of each of the two components of the equilibrium condi- 
tion are as follows: 

(a) Phase-Locking: Reduction in the phase-speed of each wave occurs through the 
interaction mechanism - the vertical velocity field produced by the distant wave acts so 
as to diminish the phase-speed of the given wave. Furthermore, if the waves are "counter- 
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propagating" (meaning, the direction of the intrinsic phase-speed, — a>i sin (A4>u) /a, is 
opposite to the background flow) , the background flow causes an additional reduction in 
the phase-speed. Both wave interaction and counter-propagation work in tandem until 
the two waves get "phase-locked" , i.e. they are stationary relative to each other. In other 
words, this means d^/dt = 0. 

(b) Mutual Growth: The phase shift at the phase-locked state, 8 nm , is a unique angle 
producing the resonant configuration. This configuration causes the two waves to grow 
equally (i.e. 71 = 72) via interaction. Eq. ( 3.13[ ) shows that equal growth rate implies 
dR/dt = 0. Furthermore, Eqs. (|3.9[) and (3.11) imply 71 = 72 = constant, meaning that 



the wave amplitudes grow at an exponential rat^f] This exponential growth explains 
why the equilibrium point (i? nm ,0 nm ) corresponds to the growing normal-mode of the 
discrete spectrum. 

WIT shows that the individual waves grow when ^ 4> ^ tt and decay when — tt ^ 
4> 0. The largest possible growth, also known as the optimal growth, occurs when 
4> = 7r/2. These results are in accordance with the analysis of Heifetz & Methven (2005). 



Using the Generalized stability theory of Farrell & Ioannou ( 1996 ) , Heifetz &; Methven 
( |2005[ ) showed that the optimal growth in a barotropic shear layer occurs when the two 
counter-propagating Rossby waves are in quadrature. 

WIT also reveals that the left moving top wave and the right moving bottom wave 
eventually phase- lock (which then leads to mutual growth), provided the condition in Eq. 
(3.18) is satisfied. Any arbitrary initial condition (say R = Rq and $ = $0) finally leads to 
phase-locking, as evident from Fig. [2] As mentioned already, Eqs. (3.9 1-( 3. 14 ) describe the 
non-modal instability process. Non-modal instability signifies non-orthogonal interaction 
between the two wave modes, and is the entire process occurring prior to the phase- 
locking event. Phase-locking is the final or steady state configuration, and corresponds 
to the the growing normal-mode (R„ m ,9„ m ) of the discrete spectrum. The fact that 9 nm 
signifies the growing normal- mode configuration implies ^ 6 nm ^ tt. 

A misconception might arise for phase shifts in the range — tt ^ $ 0. In this case, 
the reader might form an impression that the instability will not appear (because both 
the waves are decaying according to Eqs. (3.9) and (3.11)). However, the wave decay 
process is temporary. The two waves continuously adjust $ so as to enter the growing 
zone ^ $ ^ tt. After reaching this zone, the two waves still continue to adjust $ until 
the steady-state value (i.e. the resonant configuration), 9 nm , is reached. This fact can be 
better understood by considering the case <&o = — tt/2 in Fig. [2|a). Although initially 
— tt ^ $ ^ 0, the value of $ eventually enters the growing range and finally attains the 
steady-state value. 



4. Homogeneous and stratified shear instabilities 

4.1. The Kelvin Helmholtz Instability 
Let us consider a piecewise linear velocity profile 

{Ui z ^ zi 

Sz z 2 s? z sC zi 
U 2 z < z 2 



(4.1) 



f There are systems where phase-locking does not produce exponential growth. For example, 
stable barotropic and baroclinic modes result from t he phase-locking betwee n deep water surface 
gravity and internal gravity waves; see Chapter 7 of |Kundu fc Cohen (20041, Pg. 259-261. Shear 
is always absent in such systems. 
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0.25 



7 




(b) 



Figure 3: (a) The setting leading to the Kelvin-Hclmholtz instability. The velocity profile 
in Eq. ( |4.2[ ) is shown on the left, while the vorticity waves (marked by "V") are shown 
on the right, (b) Linear stability diagram of the Kelvin-Hclmholtz instability (7 denotes 
the modal growth rate). 



This profile is a prototype of barotropic shear layers occurring in many geophysical and 
astrophysical flows (Guha et al. 2013[ ). It supports two vorticity waves, one at z\ and 



the other at z 2 . The shear S = (Ui — U2)/{z\ — z 2 ). We nondimensionalize the problem 
by choosing a length scale h = {z\ — z 2 )/2 and a velocity scale AU = (Ui — U2)/2. In 
a reference frame moving with the mean flow U = (Ui + U2)/2, the non-dimensional 
velocity profile becomes 

{1 Ol 
z -HzO (4.2) 
-1 zsC-l 

where both U and z are now non-dimensional quantities. This profile, along with the 
vorticity waves, is shown in Fig.[3|a). The top wave is left moving while the bottom wave 
is right moving. Both the waves counter-propagate, i.e. move in a direction opposite to 
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the background flow. The wave interaction and subsequent instability mechanism can be 
understood from WIT, see <|3} 



The classical normal-mode based linear stability analysis of the profile in Eq. (4.2) was 



first performed by Rayleigh ( 1880 1. He showed that if the non-dimensional wavenumber a 
is in the range < a ^ 0.64, the flow is unstable; see Fig.|3jb). Thus, the piecewise linear 
profile and the ensuing instability are often referred to as the "Rayleigh's shear layer" 
and "Rayleigh's shear instability" respectively. However we will address the latter as the 
"Kelvin-Helmholtz instability (KH)", following the wider acceptance of this terminology 

The non-modal analysis 



in the stratified shear layer community (Carpenter et al. 



of the piecewise linear profile was performed in detail by 



Methven (2005). Following the footsteps of Bretherton (1966) and Hoskins et al. (1985), 



2012) 



Heifetz et al. (1999); Heifetz & 



Heifetz and co-authors were able to put forward a comprehensive mechanistic picture of 
KH (in rotating frame) in terms of counter-propagating Rossby wave interactions. By 
using the Generalized Stability Theory (Farrell & Ioannou 19961, Heifetz & Methven 



(2005) showed how wave interaction leads to optimal growth in shear layers. 



Here we study the KH problem in terms of WIT, i.e. Eqs. ( 3.9 )-( 3.12 1. Since the two 
waves involved in the KH problem are vorticity waves, we substitute Eq. (2.14) in the 
WIT equation-set and after non-dimensionalization we obtain 



7i 



72 





sin $ 


2R 




1 


1 






~~ 2a 




-e~ 2a sin$ 


2 





-2d 



COS $ 



ca = -1 + — [1 - Re- 2a cos $1 
2a J 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



Eqs. (4.3I-&6J) are isomorphic to Eqs. (14a)-(14d) of Heifetz et al. (1999) and homo- 



morphie to Eqs. (7a)-(7d) of Davies & Bishop (1994). These two referenced equation-sets 



describe edge wave interactions in two different types of rotating physical systems. While 



the one described by Heifetz et al. ( 1999 ) shows how CRW interactions lead to barotropic 



shear instability, the equation-set formulated by Davies & Bishop (1994]) shows how baro- 
clinic instability is produced through the interaction of temperature edge waves of the 
Eady model. Furthermore, Heifetz et al. fl999) showed that their set of equations is 

p94|) . 



homomorphic to that of Davies & Bishop 



Eqs. (4.3 )-(4.6 ) demonstrate how wave interaction causes amplitude growth and phase- 
speed modification of the individual vorticity waves, thereby leading to KH. The fact that 
the wave interaction modifies the phase-speed of a vorticity wave can be understood by 



comparing Eq. (4.4) and Eq. (4.6) with the non-dimensional form of Eq. (2.16) (nondi- 
mensionalization means substituting 5 = 1 and Ui = 1 or —1 in Eq. (2.16)). 



The generalized non-linear dynamical system given by Eqs. (3.13)-(3.14) in this case 
translates to 



dR _ 1 , 



R 2 ) e~ 2a sin$ 



~dt 



= (2a-l) 



R 



e 2a cos <I> 



(4.7) 
(4.8) 
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Figure 4: Phase portrait of Kelvin-Hclmholtz instability corresponding to a = 0.4. The 
system has two equilibrium points - one unstable (o) and the other stable (•). $ is the 
phase difference between the lower and upper waves, while R represents the ratio of the 
upper wave amplitude to the lower wave amplitude. 



The equilibrium points of this system are (R nm , ±^m)i where 



Rn 



O^^cos- 1 [(l-2a)e 2Q ] 



(4.9) 
(4.10) 



The phase portrait is shown in Fig. [4j It confirms that the the dynamical system is 
indeed of source-sink type, as predicted in <|3j 



The necessary and sufficient condition for instability expressed via Eq. (3.181 in this 
case translates to 



- 1 < (1 - 2a) e 2a ^ 1 implying O^a^ 0.64 



(4.11) 



The range of unstable wavenumbers obtained from the above equation corroborates 
Rayleigh's normal-mode analysis. 

Raylcigh also found the wavenumber of maximum growth to be a max — 0.4. This value 
can be verified through WIT by imposing the normal-mode condition and maximizing 
7x or 72 with respect to a. 

The fact that KH develops into a standing wave instability can be verified by applying 



the normal- mode condition in Eqs. (4.4) and (4.6). Performing the necessary steps we 
find c\ — C2 — 0, i.e. the waves have become stationary after phase-locking. In this 
configuration, the waves start to grow exponentially. Hence the shear layer grows in size. 
The growth process eventually becomes non-linear, and the shear layer modifies into 
elliptical patches of constant vorticity ( Guha et aZ.|2013 |. 
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U(z) p(z) 




(a) 



J 




(b) 

Figure 5: (a) The setting leading to the Taylor instability. The velocity and density 
profiles in Eq. (4.13) are shown on the left, while the gravity waves "G" are shown on 
the right, (b) Linear stability diagram of the Taylor instability. The contours represent 
the growth rate. 



4.2. The Taylor Instability 
Let us consider a uniform shear layer with two density interfaces 



U(z) = Sz and p (z) 




(4.12) 



The shear S is constant. We choose Ap/2 as the density scale, h = (zi — Z2)/2 as the 
length scale, and thereby nondimensionalize Eq. (4.12 1. The physical state of the system is 
determined by the competition between the density stratification and the shear, the non- 
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dimensional measure of which is given by the Bulk Richardson number J = g /(hS 2 ), 
where g = g(Ap/2)/ p is the reduced gravity, and po is the reference density. The 
dimensionless velocity and density profiles therefore become 



U(z) = z and p (z) = 




(4.13) 



This flow configuration is shown in Fig. [5|a). Contrary to the conventional notion that 
gravitationally stable density stratified flows are usually stable, Taylor ( 1931 ) put forward 



the flow given by Eq. (4.13) and showed it to be linearly unstable. The interplay between 



the background shear and the gravity waves existing at the density interfaces produce 
the destabilizing effect. The ensuing instability is thus known as the "Taylor instability". 
Taylor found that for each value of J, there exists a band of unstable wavenumbers 
(and vice-versa), shown in Fig. ^b). This unstable range is given by (see Eq. (2.154) of 
|Sutherland| p0l5| ) 

< J < ^ ^ ( 4 - 14 ) 



1 



1 



Caulfield (1994), and more recently Carpenter et al. (2012), have described Taylor 



instability in terms of wave interactions. As discussed in {2.2 each density interface 
(located at z = 1 and z = — 1) supports two gravity waves. The interaction between the 
left moving gravity wave at the upper interface and right moving gravity wave at the 
lower interface leads to Taylor instability. 



To understand Taylor instability in terms of WIT, we substitute Eq. (2.20) in Eqs 



( 3.9 )-( 3.12 ). After performing non-dimensionalization, we obtain 



7i 



J 



2R(l + c 2 ) 



e 2a sin $ 




J ( , R 9 

— 1 e~ 2a cos<I> 

2a V P 



(4.15) 

(4.16) 
(4.17) 

(4.18) 



Here j3 = CU2/0J1 = (1 — ci)/(l + C2) is the frequency ratio, and is therefore a positive 



quantity. From Eqs. (4.16) and (4.18) we construct a quadratic equation for f3: 



f3 2 + f5e- 2a cos$ 



- R - 1 =0 



(4.19) 



Only the positive one is relevant. 

The coupled nature of Eqs. (4.16) and (4.18) makes it more complicated than the KH 
problem. The non-linear dynamical system in this case is given by 
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Figure 6: Phase portrait of Taylor instability corresponding to an unstable combination 
of a and J. Here a = 0.2 and J = 0.7264. 



ip. Substituting this value in Eq. (4.191 gives /3 = 1 



At phase-locking R = R nm 
Therefore R nm = 1 and c\ = c-i — at resonance. This implies that Taylor instability, 
like KH, also evolves into a standing wave instability. Although this fact is previously 
known, WIT demonstrates why this is the case. The non-linear manifestation of Taylor 



instability is quite similar to KH. Caulfield et al. ( 1995 ) has experimentally shown that 



this instability evolves into billowing structures. 
The phase-shift 8 nm is evaluated from Eq. (4.21): 



= cos 



The necessary and sufficient condition for Taylor instability is given by 

2a , 2a 



2a\ 

— \e 2a ^l implying - - ^ 



1 



-2a 



(4.22) 



(4.23) 



The latter result corroborates the classical normal-mode result given in Eq. (4.14|. 



4.3. The Holmboe Instability 
Let us consider the following velocity and density profiles 



U(z) 




and p (z) 




(4.24) 




J 




cx 

(b) 

Figure 7: (a) The setting leading to the Holmboe instability. The velocity and density 
profiles in Eq. (4.25) are shown on the left, while the vorticity wave "V" and the gravity 



wave "G" are shown on the right, (b) Linear stability diagram of the Holmboe instability. 
The contours represent the growth rate. 



We nondimensionalize Eq. (4.24) exactly like the Taylor problem, which gives us the 



dimensionless velocity and density profiles: 



U(z) = 




JO z ^0 

^ P[Z)= 2 z^O 



(4.25) 



The vorticity interface at the top supports a vorticity wave, while the density interface at 
the bottom supports two gravity waves. The interaction between the left moving vorticity 
wave at the upper interface and the right moving gravity wave at the lower interface leads 
to an instability mechanism, known as the "Holmboe instability" . The corresponding flow 
setting is shown in Fig. . 
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Figure 8: Phase portrait of Holmboc instability corresponding to a = 1 and J = 0.5. 



Holmboe ( 1962 ) was the first to consider the instability mechanism resulting from the 



interaction between vorticity and gravity waves. In his actual problem, Holmboe consid- 



ered a flow setting more complicated than Eq. (4.25). His problem consisted of a velocity 
profile given by Eq. (4.2), however the density profile is the same as that in Eq. (4.25). 



Holmboe performed a linear stability analysis and showed that in addition to the conven- 
tional KH mode, there is another mode of instability - the Holmboe mode. Unlike the KH 
mode, the Holmboe mode is characterized by traveling waves. Presence of this unstable 
mode reveals that stable density stratification can also have a destabilizing influence. 
This aspect of Holmboe instability is much like the Taylor instability. Recent non-modal 



analysis by Constantinou & Ioannou (20111 has shown that Holmboe instability is sus- 



ceptible to substantial transient growths. Such growths especially occur for parameter 
values for which there is no instability but are close to the stability boundary. 

Analyzing the "authentic" Holmboe instability in terms of WIT implies considering the 
interaction of three waves - two vorticity waves and a gravity wave. An extended version 
of WIT can handle this problem, however this will not be considered in this papei|f] 
Baines & Mitsudera ( |1994[ ) simplified Holmboe's problem by introducing the profile in 
Eq. (4.25). This allows studying the interaction of a vorticity and a gravity wave, and 
is therefore suitable for this paper. Linear stability analysis shows that corresponding to 
each value of J, there exists a band of unstable wavenumbers. This is shown in Fig. [Tj^b) . 
The stability boundary has been evaluated in Appendix [A) 

In order to understand Holmboe instability in terms of WIT, we substitute Eq. ( |2.14| 



f Stratifie d shear layer inst abilities resulting from the interaction of multiple waves have been 
addressed by Caulfield 1 1994 k however he limited his study to the normal-mode waveform. 
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and Eq. (|2.20|) in Eqs. (|3.9|)-(|3.12|). After performing non-dimensionalization, we obtain 

(4.26) 



7i 



J 



1 ( 1 

a \ z 



J 



e Q cos $ 



1? 



72 



C'2 



: sin $ 



(-i?e" Q cos $ + ^R 2 e- 2a cos 2 $ + 16a j) 



(4.27) 

(4.28) 
(4.29) 



Like the Taylor case, this equation-set is also of coupled type. The non-linear dynamical 
system in this case is given by: 

4a J R 2 " 
-Re~ a cos $ + y / R 2 e~ 2a cos 2 $ + 16a J 2 



dR 
~dt 
d<$> 
~dt 



' sin $ 



(4.30) 



1 



a 



(1 - Re- a cos$) + 



AaJ 



-Re~ a cos $ + \ / R 2 e- 2a cos 2 $ + 16a J 



e~ Q cos$ 



R 



(4.31) 



The equilibrium points of this system are (i? nm , =b# nm ), where 

Rnm — 



\ 



1 - 2a + \ 32aJ + (1 - 2a) 



i? 



2 

nra 



2 

1 - 2a 



2il„ 



The necessary and sufficient condition for Holmboe instability is found to be 

'R* + 1 - 2a N 



1 «J 



e a ^ 1 



2Rnm 

This provides the range of J leading to Holmboe instability, and is as follows: 



1 

2A 



-B - VB 2 - 4AC 1 ) < J < (-B + VB 2 - 4AC 



(4.32) 
(4.33) 

(4.34) 
(4.35) 



where 



A = 16a 2 
B = 



8 (2a- If + 36 (2a- 1) e 

-2a 



-2a 



Tie 



-4a 



C = (2a - 1 + er 2a ) (2a - l) 3 

Eq. (4.35) corroborates the normal- mode result given in Appendix |A| 

The phase portrait of Holmboe instability, corresponding to an unstable combination 
of a and J, is shown in Fig. [HJ This phase portrait is slightly different from Taylor 
and KH, because R nm ^ 1 in this case. Another feature of Holmboe instability is that, 
unlike Taylor and KH cases, its phase-speed is non-zero at the equilibrium condition. 
This phase-speed is found to be 

2 J 4 J 

ci = c 2 - — = , (4.36) 



l-2a+ J32aJ+ (1 - 2a)' 
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In the limit of large a and J, the two phase-locked waves move with unit speed to the 
right. 

5. Summary and Conclusion 

Shear instability plays a crucial role in atmospheric and oceanic flows. In the last 50 
years, significant efforts have been made to develop a mechanistic understanding of shear 
instabilities. Using idealized velocity and density profiles, researchers have hypothesized 
that the resonant interaction between two counter-propagating linear interfacial waves is 
the root cause behind exponentially growing instabilities in homogeneous and stratified 
shear layers. Support for this claim has been provided by considering interacting vorticity 
and gravity waves of the normal-mode form. 

This paper is devoted to investigating the wave interaction problem in a generalized 
sense. The governing equations of hydrodynamic instability in homogeneous and stratified 
shear layers have been derived without imposing the wave type, or the normal-mode 
waveform. We refer to this equation-set as the Wave Interaction Theory (WIT). Using 
WIT we have shown in Fig. [2] that two counter-propagating linear interfacial waves, 
having arbitrary initial amplitudes and phases, eventually resonate, provided they satisfy 
the necessary and sufficient condition for resonant interaction (Eq. ( |3.18[ )). Resonance 
makes the interfacial waves, and therefore the shear layer, to grow at an exponential rate 
- leading to the normal-mode instability. The condition for resonant interaction provides 
the range of unstable wavenumbers causing exponential growth. By considering three 
different types of shear instabilities - Kelvin Helmholtz, Taylor and Holmboe, we have 
shown that the resonant condition in each case matches the predictions of the canonical 
normal-mode based linear stability theory. 

In addition to perceiving the shear instability problem in terms of two interacting linear 
waves, WIT also provides an alternate perspective - understanding shear instability in 
terms of dynamical systems. According to WIT, an unstable shear layer represents a 
non-linear dynamical system of the source-sink type, source and sink being the two 
equilibrium points. The equilibrium condition of the dynamical system is analogous to 
the resonant condition of the two-interacting-waves system. In terms of canonical linear 
stability theory, the source and the sink respectively correspond to the decaying and the 
growing normal-modes of the discrete spectrum. 

The most important aspect of WIT is that it provides a non-modal description of 
idealized shear instabilities. Non-modal instability signifies non-orthogonal interaction 
between the two wave modes, and is the entire process occurring prior to resonance. 
Rapid transient growth is a key feature of non-modal instability process. WIT shows 
that optimal growth occurs when the two waves are in quadrature. The phenomenon of 
transient growth has been studied very briefly in this paper; thorough research is required 
in future for its detailed understanding. 

Although we have limited our study to the interaction between two waves, the Taylor 
and Holmboe profiles actually involve multiple wave interactions, which we have ne- 
glected. Each density interface in these profiles supports two gravity waves, out of which 
only the counter-propagating wave has been considered. Noting that two wave interac- 
tions are sufficient to produce the normal mode characteristics of Holmboe and Taylor 
instabilities, it can be argued that the inclusion of co-propagating gravity wave would 
have been unnecessary. However this wave might have some effect during the initial 
interaction (non-modal) stages, which needs to be studied in future. 

Finally we focus on the implications of using idealized velocity and density profiles. 
These idealizations have allowed us to concentrate only on the discrete spectrum dynam- 



22 



A. Guha and G.A. Lawrence 



ics and understand hydrodynamic instability in terms of interfacial wave interactions. 
Real profiles are always continuous, which indicates the importance of the continuous 
spectrum. Using Green function technique, Heifetz & Methven (2005) has shown that 
the continuous spectrum dynamics in a smooth, homogeneous shear layer can be under- 



stood in terms of infinite number of interacting vorticity (Rossby edge) waves. Harnik 



et al. ( 2008 ) used the same approach to understand the normal- mode continuous spec- 



trum of smooth, stratified shear layers. Summarizing, our two-wave-interaction problem 
gets modified into an infinite-wave-interaction problem in the case of continuous profiles, 
however the underlying mechanism remains the same. 



Appendix A. Normal mode form of Holmboe instability 

Both interfaces in the Holmboe profile (Eq. ( |4.25 1) individually satisfy the kinematic 



condition: 



dm 

dt 

dm 

dt 



dx 
d_ 
dx 



1 - 2a 



2a 



(Al) 
(A 2) 



where tp2 is the stream function perturbation at the lower interface. This interface being 
a density interface also satisfies the dynamic condition: 



dip 2 
dx 



J dr] 2 
a dx 



(A3) 



We assume the perturbations to be of normal-mode form: rji — f]\e %ai - x ct \ m = 
^ e ia{x-ct) ^ anc j ^ 2 _ ^ 2 gia(x-ct)_ jj ere ^ ne wave speed c is generally complex. Defining 

? = [ "02 f\i 771 ] , we obtain the following eigenvalue problem: 

(M + c/)? = (A 4) 

where 

J/a 



M 





e~ a /{2a) 
(I - 2a) /(2a) 



Eq. (A4) generates the following characteristic polynomial 

2a 



1 



2a 



J -c-J 
a 



1 - 2a 
2a 2 



-2a 



J 



2a 2 







(A 5) 



(A 6) 



This equation produces complex conjugate roots only when the discriminant is negative. 
Since the presence of complex roots signify normal-mode instability, negative values of 
the discriminant is of our interest. The discriminant (D) in this case is given by: 



D = 16a 2 J 2 -aJ 



l(2a- 1) 



36e' 



(2a ~1) + 27e~ 



(1 - 2aY (2a - 1 



Imposing the condition D < 0, we find 



I 

2A 



-B - VB 2 - 4AC) J sc i (-B + VB 2 - 4AC 



(A 7) 
(A 8) 
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where 

A = 16a 2 



B = -a 



8 (2a - l) 2 + 36 (2a - 1) e~ 2Q + 27e 



—4a 



C = (2a- 1 + e~ 2a ) (2a- 1) 



Thus Holmboe instability occurs only when the condition in Eq. ( A 8 1 is satisfied 



REFERENCES 

Baines, P.G. k, Mitsudera, H. 1994 On the mechanism of shear flow instabilities. J. Fluid 
Mech. 276, 327-342. 

Bretherton, F. P. 1966 Baroclinic instability and the short wavelength cut-off in terms of 

potential vorticity. Q. J. Roy. Meteor. Soc. 92 (393), 335-345. 
Cairns, R.A. 1979 The role of negative energy waves in some instabilities of parallel flows. J. 

Fluid Mech. 92, 1-14. 

Carpenter, Jeffrey R., Tedford, Edmund W., Heifetz, Eyal & Lawrence, Gre- 
gory A. 2012 Instability in stratified shear flow: Review of a physical interpretation based 
on interacting waves. Applied Mechanics Reviews 64 (6), 061001. 

Caulfield, CP. 1994 Multiple linear instability of layered stratified shear flow. J. Fluid Mech. 
258, 255-285. 

Caulfield, CP., Peltier, W.R., Yoshida, S. & Ohtani, M. 1995 An experimental inves- 
tigation of the instability of a shear flow with multilayered density stratification. Phys. 
Fluids 7, 3028-3041. 

Constantinou, Navid C. & Ioannou, Petros J. 2011 Optimal excitation of two dimensional 

holmboe instabilities. Phys. Fluids 23 (7), 074102. 
Davies, HC & Bishop, CH 1994 Eady edge waves and rapid development. J. Atmos. Sci. 

51 (13), 1930-1946. 

Drazin, P.G. & Reid, W.H. 2004 Hydrodynamic Stability, 2nd edn. Cambridge University 
Press. 

Farrell, B. 1984 Modal and non-modal baroclinic waves. J. Atmos. Sci. 41 (4), 668-673. 

Farrell, B.F. & Ioannou, P.J. 1996 Generalized stability theory, part i: Autonomous oper- 
ators. J. Atmos. Sci. 53 (14), 2025-2040. 

Goldstein, S. 1931 On the stability of superposed streams of fluids of different densities. Proc. 
R. Soc. Lond. A 132, 524-548. 

Guha, Anirban, Rahmani, Mona & Lawrence, Gregory A. 2013 Evolution of a barotropic 
shear layer into elliptical vortices. Phys. Rev. E 87, 013020. 

Harnik, N., Heifetz, E., Umurhan, OM & Lott, F. 2008 A buoyancy- vorticity wave in- 
teraction approach to stratified shear flow. Journal of the Atmospheric Sciences 65 (8), 
2615-2630. 

Heifetz, E., Bishop, CH & Alpert, P. 1999 Counter-propagating Rossby waves in the 
barotropic Rayleigh model of shear instability. Q. J. R. Meteorol. Soc. 125 (560), 2835- 
2853. 

Heifetz, E., Bishop, CH, Hoskins, BJ & Methven, J. 2004 The counter-propagating rossby- 

wave perspective on baroclinic instability, i: Mathematical basis. Quarterly Journal of the 

Royal Meteorological Society 130 (596), 211-231. 
Heifetz, Eyal & Methven, John 2005 Relating optimal growth to counterpropagating Rossby 

waves in shear instability. Phys. Fluids 17 (6), 064107. 
Holmboe, J. 1962 On the behavior of symmetric waves in stratified shear layers. Geofys. Publ. 

24, 67-112. 

Hoskins, B.J., McIntyre, M.E. & Robertson, A.W. 1985 On the use and significance of 

isentropic potential vorticity maps. Q. J. Roy. Meteor. Soc. Ill (470), 877-946. 
Kundu, PK & Cohen, IM 2004 Fluid Mechanics. Elsevier, Boston. 

Lindzen, R.S. 1988 Instability of plane parallel shear flow (toward a mechanistic picture of how 
it works). Pure Appl. Geophys. 126 (1), 103-121. 



24 



A. Guha and G.A. Lawrence 



Orr, W. M. F. 1907 Stability or instability of the steady motions of a perfect liquid and of a 

viscous liquid. Proc. Roy. Irish Acad. A (27), 9-138. 
Rayleigh, J.W.S. 1880 On the stability, or instability, of certain fluid motions. Proc. Lond. 

Math. Soc. 12, 57-70. 

Schmid, P.J. & Henningson, D.S. 2001 Stability and transition in shear flows, , vol. 142. 
Springer Verlag. 

Sutherland, B.R. 2010 Internal gravity waves. Cambridge University Press. 

Taylor, G.I. 1931 Effect of variation in density on the stability of superposed streams of fluid. 

Proc. R. Soc. Lond. A 132, 499-523. 
Tedford, E.W., Pieters, R. & Lawrence, GA 2009 Symmetric holmboe instabilities in a 

laboratory exchange flow. J. Fluid Mech. 636, 137-153. 
Thorpe, S.A. 1973 Experiments on instability and turbulence in a stratified shear flow. J. Fluid 

Mech. 61, 731-751. 

Trefethen, L.N., Trefethen, A.E., Reddy, S.C. & Driscoll, T.A. 1993 Hydrodynamic 

stability without eigenvalues. Science 261, 578-584. 
Turner, J.S. 1973 Buoyancy Effects in Fluids, first paperback edn. Cambridge University 

Press. 



